Pair production beyond the Schwinger formula in time-dependent electric fields 
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We investigate electron-positron pair production in pulse-shaped electric background fields using 
a non-Markovian quantum kinetic equation. We identify a pulse length range for subcritical fields 
still in the nonperturbative regime where the number of produced pairs significantly exceeds that 
of a naive expectation based on the Schwinger formula. Prom a conceptual viewpoint, we find 
a remarkable quantitative agreement between the (real time) quantum kinetic approach and the 
(imaginary time) effective action approach. 

PACS numbers: 12.20.Ds, ll.15.Tk, 



I. INTRODUCTION 

It is a long-standing prediction of quantum electrody- 
namics (QED) that its vacuum is unstable in the presence 
of strong electric fields and decays by emitting electron- 
positron pairs [l|, y • The decay of the QED vacuum - the 
so-called Schwinger mechanism - has often been analyzed 
within equilibrium quantum field theory, even though it 
is a far-from-equilibrium, time-dependent phenomenon. 
For resolving the dynamics of the production process, ki- 
netic theory including a source term for electron-positron 
pair production provides an appropriate approach. 

According to the nonperturbative Schwinger formula 
[3|, which is strictly valid only for spatially homoge- 
neous and static electric fields, a sizeable production 
rate requires field strengths of the order Ecr = 1.3 • 
10^® V/m. The study of corrections due to spatial or 
temporal inhomogeneities has a long history which in- 
cludes exact solutions J, 01 , various semiclassical meth- 
ods [0, 0, H S E, [ll|, and functional techniques [l2| 
as well as Monte Carlo simulations [13|. Within ki- 
netic theory, first approaches used an instantaneous phe- 
nomenologically motivated source term based on the 
Schwinger formula, which served as a model for both 
quark-antiquark pair production in chromoelectric fiux 
tubes and electron-positron pair production even in the 
case of time-dependent electric fields |14. iisl. UM. NLTI. [ia|. 
On the other hand, for spatially homogeneous and time- 
dependent electric fields, a rigorous connection between 
kinetic theory and a mean-field approximation to QED 
has been estabhshed [13, [20, [2l|. It was shown that the 
true source term for electron-positron pair production 
has intrinsically non-Markovian character. 

In recent years, special interest has been laid on the 
investigation of electron-positron pair production in al- 
ternatin g [ 22. [23 . [24 | and pulse-shaped electric fields 
|25l . [20 . [271 1 within the kinetic approach including the 
non-Markovian source term, with the field-current feed- 
back due to Maxwell's equation taken into account as 
well. Because of the rapid development of laser technol- 
ogy during the past years, physicists pin their hopes on x- 



ray free electron laser systems and optical high-intensity 
laser systems to experimentally observe the Schwinger 
mechanism within the next decades. 

In this investigation, we explore the connection bet- 
ween the non-Markovian and the Schwinger-like source 
term for a pulse-shaped electric field in order to look for 
a new pair-production behavior which goes beyond the 
instantaneous Schwinger approximation. 



II. PARTICLE CREATION IN ELECTRIC 
BACKGROUND FIELDS 

We consider a spatially homogeneous, time-dependent 
AbeHan vector potential in the temporal gauge Af^{t) ~ 
{0,A{t)), with its spatial part defining the z-direction: 
A{t) ~ {0,0, A{t)). The corresponding electric back- 
ground field points into the z direction as well. 



m 



dA{t) 
dt 



(0,0, E{t)) 



(1) 



whereas the corresponding magnetic background field 
vanishes: B(t) = V x A{t) ~ 0. In our calculations, 
we choose an electric background field of the form 



E{t) = 



En 



cosh^ (t/r) 



(2) 



The field reaches a maximum magnitude of Eq a.t t = 
and switches on and off exponentially near t « ±4t, 
where the electric background field drops to a value of 
approximately one per mille of the maximum magnitude. 
The temporal-gauge vector potential A{t) giving rise to 
this electric background field is given by 



A{t) = -EoT [1 + tanh (t/r)] 



A. Quantum kinetic equation 



(3) 



A key quantity in the description of the pair- 
production process in the electric background field is the 



single-particle momentum distribution function f{k,t). 
It can be computed from a quantum Vlasov equation 
with a source term for electron-positron pair production. 
It receives a physical meaning as distribution function of 
real particles as we go to asymptotic times t ^ ±00 when 
the electric background field vanishes. In our opinion, a 
similar interpretation at intermediate times [2J| is diffi- 
cult to justify. In addition to the source term, we would in 
general have to take collisions and the field-current feed- 
back due to Maxwell's equation into account; however, in 
the field strength regime we are mainly interested in, i.e. 
the subcritical field strength regime with i?o ~ O.lE'cr, 
these effects can be neglected [25|, and we obtain 



dr^(fc,t) 

di 



1 eE{t)ex_ 



dt 



, eE{t')e^ 

LU^{k,t') 



(4) 



i-2r^{k,t') 



cos 



2 / dTu;{k,T) 
t' 



with e the electric charge, k = {k±,k3) denotes the 
canonical three-momentum vector with p|| (t) = k^ — 
eAit) being the kinetic momentum in the z direction, 
e\ = m? + fc^ is the transverse energy squared and 
u}'^{k,t) = e']_ + p||(t)^ characterizes the total energy 
squared. This equation, which is valid for both helici- 
ties, has a non-Markovian character for two reasons: the 
factor [1 — 2/(fc, t)] arising from quantum statistics takes 
the time history of the distribution function into account, 
whereas the cosine factor accounts for the time history 
of the electric background field. In order to simplify the 
numerical treatment of this equation, we may introduce 
the auxiliary quantities 



iM 



{k,t) = 



„nM 



(fc,<) 



00 W2(fc,t') 

-2r^\k,t') 

00 W2(fc,t') 



dr Lo{k, t) 



dr uj{k, t) 



(5) 



(6) 



and reformulate the integro-differential equation (JU as a 
coupled system of first order differential equations: 



dt 

d 

dt' 

d 

dt' 



r\k,t) 



1 eE{t)e±_ 
^u^{k,t) 



,nM 



{k, t) , 



,nM 



,nM 



{k,t) 
{k,t) 



2uj{k,ty'^{k,t) 

eE{t)e± 



l-2r''\k,t) 



u^k.t) 
-2cj{k,t)u''^{k,t) 



(7) 

(8) 
(9) 



While the non-Markovian character of Eq. (j4]) due to 
the cosine factor is crucial for the electron-positron pair 



creation process, we anticipate that memory effects due 
to the statistical factor become only important for very 
strong electric background fields of the order of Eq « 
Ecr- Then, the Compton time t^ = l/m is of the order 
of the time scale of electron-positron pair creation [2a|, 
and a sizeable distribution localized in phase space is 
accumulated such that statistical effects such as Pauli 
blocking become relevant. However, in the field strength 
regime we are interested in, Eq « O.li^cr, neglecting the 
statistical factor for f{k, t) ^ 1 yields the low-density 
approximation: 



s^"A" 



1 e£:(t)ej 

^u;^k,t) J -00"^ u;^k,t') 



' dt'^^m^ (10) 



X COS 



dTUj{k, t) 



This equation can be solved either by direct integration 



.r(fc,i) 



^,,eEit")e, 



t" 



LU^k,t") 

, ,eE{t')e±_ 

dt' V^-^cos 

uj^{k,t') 



(11) 



t" 



dTUj{k, t) 



or, as before, by reformulating it as a coupled system of 
first-order differential equations: 



^^u"'{k,t) = 2cj{k,ty^{k,t) , 



(12) 

(13) 



dt 



w'^(fc,t) 



eE{t)ej_ 

UJ^ik,t) 



2u{k,t)v}'^{k,t) . (14) 



An important observable is given by the asymptotic 
particle number density n[e^e~] which is obtained as an 
integral of the asymptotic distribution function f{k, cxd) 



i[e+e-] = 2 



d^k 

(2^ 



-/(fc,oo) , 



(15) 



with the factor of 2 arising from spin degeneracy. 



B. Effective action approach 

We briefiy summarize the main results obtained within 
the effective action approach to electron-positron pair 
production for later use. Within this approach the 
vacuum-to- vacuum persistence amplitude in the presence 
of an external electromagnetic field A^^ is represented as 



(0|0)^-CXp[iS'eff] 



(16) 



with 5'oft denoting the effective action. At weak fields, 
the imaginary part of the effective action is an estimate 



of the number of produced electron-positron pairs per 
unit volume ^,^, 



rf^\e+e 



2Im[5o 



(17) 



For the present electric field Eq. |[2|, an exact integral 
representation of the imaginary part of the effective ac- 
tion can be found 0, @| : 



Im[5eft] = -gi^ y d^fc In [(1 - e-^"+)(l - e-'^"-)] 



(18) 



with Vl± = '''{a^ + a_ ± 2ei?oT) and a± 



C. Instantaneous approximation 

Results concerning electron-positron pair production 
in time-dependent electric fields E{t) within both the 
quantum kinetic approach and the effective action ap- 
proach were derived long after Schwinger had found an 
analytic expression for the vacuum decay rate in a spa- 
tially homogeneous but static electric field Eq pointing 
into the z direction: 



Im [Scs] 



'E-^ 



87r3 



V- 

"^-^ n 

n=l 



1 



exp 



eEn 



(19) 



As an ansatz, this result was assumed to approximately 
hold for arbitrary time-dependent electric background 
field E{t) as well. This served as a starting point for 
formulating a Vlasov equation for the single-particle mo- 
mentum distribution function including an instantaneous 
phenomenologically motivated source term for electron- 
positron pair production |14. fLTI. [la| : 



di 



fik.t) 



X In 



eE{t)S{k3 -eA{t)) 



1 — cxp — - 



eE{t) 



(20) 



The instantaneous source term is chosen in such a way 
that Eq. ifTO]) is recovered [with Ei^ being replaced by 
E{t)] upon integrating over the whole momentum space 
and by including a factor of 2 for spin degeneracy. Of 
course, using the instantaneous source term is not jus- 
tified in the presence of a rapidly varying electric back- 
ground; deviations of the full result from this approxima- 
tion for time-dependent fields E{t) point to new proper- 
ties of pair production beyond the Schwinger formula, 
which will be investigated in the next section. 



III. NUMERICAL RESULTS 

In our studies, we compare n[e+e~] obtained by means 
of the numerical solution of both the non-Markovian 
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FIG. 1: Logarithmic plot of the number density n[e'^e~] cal- 
culated for Eq — 0.1-Ecr as a function of the temporal width r 
in units of [2 • 10~^ eV~^]. Dashed line: solution obtained for 
the non-Markovian equation ([4|) and within the effective ac- 
tion approach Eq. iflSJI . Dotted line: solution obtained within 
the low-density approximation Eq. I|10p . Solid line: solution 
obtained in the instantaneous approximation Eq. I|19p. 



equation ([4]) and the low-density approximation lfTO|) 
with the exact effective action result Eq. lfT8|) and the 
instantaneous approximation lfT9| . in order to look for 
new pair-production properties in time-dependent elec- 
tric background fields E{t) which are unexpected from 
the Schwinger formula (fT9|) . We also check the validity 
of the low-density approximation in this field strength 
regime. The numerical solutions are found by first solv- 
ing the coupled system of first-order differential equa- 
tions lUl)-® and lO-lIISl), respectively. We use an 
adaptive Runge-Kutta method and then perform the mo- 
mentum integral Eq. (fT5| on a finite momentum grid. 

Obviously, the asymptotic particle number density 
within the instantaneous approximation should be di- 
rectly proportional to the temporal width t of the electric 
pulse; cf. Eq. ([19]): 



n^\e+e-] = 



4,^3 

1 
cosh'*(i' 



dt' 



(21) 



■ cxp 



7rTO^cosh^(t') 



A. Pulse length and field strength dependence 

As a first aim, we explore the asymptotic particle num- 
ber density n[e+e~] reached in an electric background 
field with Eq = 0.1-Ecr as a function of the temporal 
width r ranging from 2 • IQ-^eV"^ to 2 • 10-^eV"\ 
cf. Fig. [l] and Table HI The shortest time scale has 
been chosen such that the Keldysh adiabaticity parame- 
ter 7 = Ecr / {EniT) ~ 1, separating the nonperturbative 
Schwinger regime 7^1 from the perturbative multipho- 
ton region 7 S> 1. 

First we note that the results obtained by means of 




FIG. 2: Momentum space plot of the asymptotic distribution 
function /""^(fc, oo) calculated for Eo = O.liScr and r — 2 ■ 
10~ eV~ ~ 10 -tc. It is an advantage of the quantum kinetic 
approach that the position and the width of the peak of the 
distribution function in momentum space can be calculated 
very precisely. 



T = 2 • 10~^ eV~ w 10 • tc, with t^ corresponding to the 
Compton time. The Schwinger result is obtained again 
for temporal widths t > 70 • ic- Our results confirm the 
observation that time variations of the field generically 
increase the pair-production yield [71, [ill . 

As a second aim, we explore the asymptotic parti- 
cle number density n[e~^e~] reached in an extreme short 
pulse, i. e. r = 2 • 10~^eV~^ « 10 ■ tc, as a function of 
the field strength of the electric background field with Eq 
ranging from O.l^cr to E„, cf. Fig. [Sand Table H The 
numerical results indicate that the instantaneous approx- 
imation significantly underestimates the reachable parti- 
cle number density for subcritical field strengths, whereas 
the Schwinger result is obtained again for field strengths 
of the order oi Eq > 0.5£'cr where the Keldysh parameter 
drops below 7 ~ 0.2. 



B. Non-Markovian solution vs low-density 
approximation 



the non-Markovian equation ^ and the effective action 
approach Eq. (fTSl) coincide up to the third digit, with 
the numerical deviation arising from the finite numer- 
ical accuracy. Hence, we have shown the equivalence 
between the (imaginary time) effective action approach 
and the (real time) quantum kinetic approach. Moreover, 
the quantum kinetic calculation contains in addition to 
the particle number density ^[e+e"] useful information 
about the momentum space distribution /(fc, 00) of the 
electron-positron pairs, cf. Fig. [2l 

Second, comparing the results obtained by means of 
the non-Markovian equation (j4]) and the instantaneous 
approximation Eq. (fT9|) indicates that the latter signifi- 
cantly underestimates the reachable particle number den- 
sity at short times. For example, deviations greater 
than a factor of 400 occur at extreme short time scales 



In Tableiniwe show the ratio n''^[e+e"]/n"'^[e+e"] in 
order to explore the validity region of the low-density ap- 
proximation. For weak and static electric fields E{t) = 
Eq < Ecr, the applicability of this approximation has 
been shown [30|. For pulse-shaped fields, the low-density 
approximation yields good results with deviations of the 
order of 10% only in the field strength regime Eq « 
O.lE'cr but becomes unreliable when reaching higher field 
strengths. This deviation increases only very slowly with 
pulse length, whereas for increasing field strengths the de- 
viation rises significantly. Although the quantitative re- 
sults within the low-density approximation get worse for 
higher field strengths, the approximation does not break 
down completely but still gives a reasonable qualitative 
picture. However, for accurate predictions of the single- 
particle momentum distribution function f{k, c») and the 
particle number density n[e"'"e~] one has to take all non- 
Markovian effects into account. Additionally, near and 



TABLE I: n"^[e+e"]/n^[e+e"] as a function of the temporal 
width r in units of [2 • 10"^ eV"^] and the ratio £ — Eo/Ea. 
Especially in the vicinity oi t ^ 10 • ic and Eq = O.l-Bcr we 
obtain a vast discrepancy. Points with deviations greater than 
10% are labeled boldface. 
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TABLE II: n''^[e+e"]/n"^[e+e"] as a function of the temporal 
width T in units of [2 ■ 10~^ eV~^] and the ratio £ = Eo/Ecr- 
The numerical results show a strong £ dependence whereas 
the influence of the pulse length on the result is very weak. 
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FIG. 3: Logarithmic plot of the number density calculated for 
T — 2 ■ 10~^ eV^^ ~ 10 ■ tc as a function of the field strength 
-Bo. Dashed line: solution obtained for the non-Markovian 
equation Q and within the effective action approach Eq. ifTSll . 
Dotted line: solution obtained within the low-density approx- 
imation Eq. lfTO|l . Solid line: solution obtained in the instan- 
taneous approximation Eq. I|20p . 



fields by means of an existing quantum kinetic formula- 
tion. A thorough comparison is made to results obtained 
from an exact effective action method and a phenomeno- 
logically motivated instantaneous approximation. 

One main result is the perfect agreement between the 
(real time) quantum kinetic approach and the (imaginary 
time) effective action approach. This is quite remarkable 
since the concepts behind these two approaches are very 
different: The quantum kinetic equation describes the 
real time nonequilibrium evolution of the quantum sys- 
tem whereas the effective action approach is an imaginary 
time equilibrium method. 

Additionally, we have shown that the instantaneous 
approximation drastically underestimates the reachable 
particle number density in the field strength regime Eq « 
O.li^cr when dealing with extreme short temporal widths 
T w 2 ■ 10~^eV~"'^ w 10 • tf.. Alongside, we have recon- 
firmed that the low-density approximation of the quan- 
tum kinetic equation gives reasonable qualitative results 
in the subcritical field strength regime. 



beyond Eq « E'er, also the field-current feedback due to 
Maxwell's equation eventually needs to be taken into ac- 
count. More detailed comparisons between different ap- 
proximations to the quantum kinetic equation including 
those with backreaction can be found in 13 ill. 



IV. CONCLUSIONS 



Our results suggest similar findings for larger pulse 
lengths but smaller field strengths, yielding similar 
Keldysh parameters 7 ~ 1. Numerical studies in this 
regime are challenging. Given the future potential of 
pulse shaping with higher harmonics, and in the light 
of suggestions for enhancing pair production with tay- 
lored pulses [32|, such investigations are highly relevant 
and clearly a worthwhile pursuit. 



We have numerically investigated electron-positron 
pair creation in short pulse-shaped electric background 
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